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We calculate the relaxational dynamical critical behavior of systems of 0(n||) © 0{n±) symmetry 
including conservation of magnetization by renormalization group (RG) theory within the minimal 
subtraction scheme in two loop order. Within the stability region of the Heisenberg fixed point 
and the biconical fixed point strong dynamical scaling holds with the asymptotic dynamical critical 
exponent 2; = — 1 where (j) is the crossover exponent and v the exponent of the correlation 

length. The critical dynamics at ny = 1 and n± = 2 is governed by a small dynamical transient 
exponent leading to nonuniversal nonasymptotic dynamical behavior. This may be seen e.g. in the 
temperature dependence of the magnetic transport coefficients. 



PACS numbers: 05.50.+q, 64.60.Ht 

I. INTRODUCTION 

In two preceding papersii^ we have considered the crit- 
ical statics and relaxational dynamics of 0(?^||) © 0{n^) 
physical systems near the multicritical point where the 
two phase transition lines of the system corresponding 
to 0(n||)-symmetry and 0(n_L)-symmetry meet. The 
space of the order parameter (OP) dimensions n\\ and 
n± decomposes in regions where the multicritical be- 
havior is described by different fixed points (FPs) - the 
0{n = ri|| + n±)- isotropic FP, the biconical FP, the de- 
coupling FP - and a region where no stable FP is found 
(run away region) (see Fig 1 in Ref.^). In the resummed 
two loop order field theoretic treatment it was found that 
for integer values of n|| and n± the biconical FP is stable 
only for a system with n|| = 1, n± = 2, and its symmetric 
counterpart. For specific initial conditions of the nonuni- 
versal parameters of the system also the 0{n = n|| + nj_)- 
isotropic FP (Heisenberg FP) might be reached. Such a 
system is physically represented by an antiferromagnet 
in an external magnetic field. The two cases mentioned 
above correspond to tetracritical and bicritical multicrit- 
ical points correspondingly. If no FP is reached the mul- 
ticritical point might be of first order, i.e. a triple point. 

The dynamics of the antiferromagnet in a magnetic 
field is quite complicated and the equations of motion 
have been formulated for the slow densities by Dohm 
and Janssen^. These Eqs. contain reversible and irre- 
versible coupling terms between the OPs (the compo- 
nents of the staggered magnetization parallel and perpen- 
dicular to the magnetic field) and one conserved density 
(the parallel component of the magnetization). In fact 
there is a second conserved density (CD) - the energy 
density - which in general should be taken into account 
but it will not be included here since in two loop order 
the specific heat exponent at the biconical FP turned out 



to be negative^ for the case n|| = 1 and n± = 2. Con- 
cerning the new static results^ a simplified dynamical 
model^ has been reconsidered^ consisting of two relax- 
ational equations for the two OP. The timescale ratio v 
between the two relaxation rates Fn and F^ introduces a 
very small dynamical transient since the dynamical FP 
lies very near to the stability boundary separating strong 
and weak dynamical scaling. The strong dynamical scal- 
ing FP is governed by v* finite and different from zero 
whereas the weak dynamical scaling FP by u"^ = or cxd 
correspondingly. 

A further step to the complete model is to include the 
diffusive dynamics of the slow CD leading to a model C 
like extension. In this extended model a new timescale 
ratio appears defined by the ratio of one of the OP- 
relaxation rate to the kinetic coefficient A of the con- 
served density to. This model has been studied in one 
loop order in RefsS.!^:^ taking into account only a part 
of dynamical two loop order terms and one loop statics. 
Here we present a complete two loop order calculations. 

The inclusion of further densities beside the OP makes 
it necessary to extend the static functional of the usual 
0'^-theory, although the OP alone would be sufficient to 
describe the static critical behavior. Such an extended 
static functional for isotropic systems (0(n) symmetry) 
with short range interaction has the form^'^ 
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Here the order parameter (j>o = 4>o{x) is assumed to be a 
n-component real vector, the symbol ■ denotes the scalar 
product. The secondary density toq = too(x) is consid- 
ered as a scalar quantity and h is the conjugated field 



2 



to TTiQ. It is chosen to have a vanishing average value 
(mo) — 0. Within statics, the above functional is equiv- 
alent to the Ginzburg-Landau-Wilson(GLW)-functional 

+ |(<?o-^o)'| , (2) 

where f is proportional to the temperature distance to 
the critical point and u is the fourth order coupling in 
which perturbation expansion is usually performed. The 
GLW-functional ^ is obtained by integrating out the 
CD, which appears only in Gaussian order in ([1]), in the 
corresponding partition function. The parameters f, u 
and 7 in ([1]) and, r and u in ([2]) are related by 



f = r + 



(3) 



The extended static functional appears in the driving 
force of the equations of motion for the OP (po and the 
CD mo. The ratio of the kinetic coefficient T in the re- 
laxation equation for the OP and the kinetic coefficient 
A in the diffusive equation for the CD defines the dynam- 
ical parameter w whose FP value governs the dynamical 
scaling of the model. 

It is worthwhile to summarize some results for model C 
at a usual critical point, where the CD can be identified 
with an energy like density and the value of specific heat 
exponent a (in any case) governs the relevance of the 
asymmetric static coupling 7 between the OP and the 
CD. Namely, this coupling is irrelevant in the RG sense 
- it vanishes at the FP - if the specific heat exponent is 
negative, i.e. the specific heat of the system does not 
diverge at the critical point. If the specific heat diverges 
then there remain two possibilities for the dynamical FP: 
either the FP value of the time scale ratio w between the 
timescale of the OP and the CD is different from zero 
and finite, or its FP value is zero or infinite. In the first 
case strong dynamical scaling with one timescale for the 
OP and the CD is realized with one dynamical scaling 
exponent z = 2 + a/v (v the exponent of the correlation 
length). In the second case weak dynamical scaling is 
present and the time scale of the OP is different from the 
timescale of the CD, both represented by a corresponding 
dynamical critical exponent. This region of the weak 
dynamical scaling FP is tiny (see e.g Fig. 1 in Ref;^). One 
should note that at the usual critical point an asymmetric 
coupling to the OP as given in Eq. ([1]) is always 'energy- 
like' independent of its physical origin. That means the 
divergence of the CD susceptibility is always described 
by the specific heat exponent a. 

In the case of a multicritical point treated here the sit- 
uation is more complicated since there are two OPs and 
a CD might couple to both of these OPs. As has been 
shown in paper I after a proper rotation (Eq. (64)) in 
the OP space temperature and magnetic like field direc- 



tions can be identified. In consequence at the multicriti- 
cal point one has to discriminate the case of energy and 
magnetization conservation. 

The paper is organized as follows. In section |TT] we ex- 
tend model C of the 0{n) symmetrical critical system to 
the case of a 0{n^\) 0{n±) symmetrical multicritical 
point as considered in papers I and II. The renormaliza- 
tion is performed in section IIIII and the field theoretic 
functions are calculated in section [TVl Then we discuss 
the possible FP and their stability in section [V] Effec- 
tive dynamical critical behavior is considered in section 
IVII followed in section IVIII by a short summary of the 
results and an outlook on further work to be done. 



II. MODEL C FOR MULTICRITICAL POINTS 

A. Static functional 

In order to describe the multicritical behavior the n- 
dimensional space of the order parameter components is 
split into two subspaces with dimensions nj_ and n|| with 
the property n^+7T,|| = n. The order parameter separates 
into 
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where (f)±o is the nj^-dimensional order parameter of the 
nj^-subspace, and (/)jjo is the n||-dimensional order param- 
eter of the n||-subspace. Introducing this separation into 
the GLW-functional ^ one obtains 
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which represents a multicritical Ginzburg-Landau- 
Wilson model. The properties of this functional con- 
cerning the renormalization, regions of stable FPs, and 
corresponding type of multicritical behavior has been ex- 
tensively discussed in paper I (see^ for earlier references) . 
The separation (j4]) has now to be performed in ([T]). The 
resulting functional is 
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+ ^i±mo(f±o ■ $±0 + ^7||"io<?||o • $\\Q - hnioj . (6) 

Integrating the contributions of the secondary density in 
the corresponding partition function, ([6]) reduces to the 
static functional (O . Relations analogous to ^ between 
the parameters of the two static functionals arise. They 
read 



^11 = ^11 ' 



u±=u±- 87^ , 



(7) 



M|| = U|| - 37|| , (8) 
Mx = Ux - 37_l7|| . (9) 



Because the partition function calculated from ^ is re- 
ducible to a partition function based on (O by integra- 
tion, the correlation functions, or vertex functions respec- 
tively, of the secondary density mo are exactly related to 
correlation functions of the order parameter. This leads 
to several relations which are important for the renor- 
malization. In particular, the average value of toq a-nd 
the two-point correlation function are defined as 



(mo) = 



1 



[mo mo) 



Bi 
1 



'D{(f)±o,<l>\\o,mo) mo e ^s-' , (10) 
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with TV^'^'' = J I?((/)xo, 0||O7 "^0) e^^^' as the normaliza- 
tion constant and 25(0X0, (/'yo, "lo) as a suitable integral 
measure. Performing the integration over mg in ()10|) and 
using Eqs.(l7])-([9]), the average value of toq reads 



(mo) = /i - 7_L ( - 
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where 0^ denotes quadratic insertions of the order pa- 
rameter. Their average values on the right hand side of 
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are now calculated with the static functional (O and 
Nsi = / 25(0X0, 0||o) e~^^'. In order to obtain (mo) = 
the conjugated external field is chosen to 



7-l( ^0io) +7|| {\<P\\o 



(14) 



Quite analogous by integrating mo in ([TT|) one obtains the 
following relation for the two-point correlation function 
of the secondary density 
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(mo mo)c = 1- 7 -r -7. 
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In (1151) where we have introduced the column matrix 



7 



7x 
711 



(16) 



The superscript indicates a transposed vector or ma- 
trix, while the subscript c on the average at the left 
hand side of ([TS]) denotes the cummulant {A B)c = 

{A B) - {A){B). The matrix 



p(0,2) f,(0,2) 

;-LJ- ^ ;X|| 
f(0,2) ^(0,2) 

^llll 



a2 1, 
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of two-point vertex functions is related to correlations of 
02-insertions. The vertex functions generally have been 
introduced in paper I (section III Renormalization) , and 
especially the matrix (jl7p (renormalized counterpart) in 
Eq.(83) therein. A third important relation can be ob- 
tained by differentiating the average value pO)l by h at 
fixed parameter l\fa,Ua,"ya- In Af„ = fa — fa^c the 
shift of the critical temperature has been taken into ac- 
count (for more details see Appendix A in paper I) . As a 
result one obtains 



d_ 

dh 



{mo{x)) 



dx' {mo{x) mo{x'))c ■ (18) 



Due to relation (fT4| the external field is function of 
Afa- The /i-derivative in p8p can be rewritten as Afa- 
derivatives. Finally one obtains 
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where we have defined 
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(19) 
(20) 



All static vertex functions, for the order parameter as 
well as for the secondary density, may be calculated with 
([S]) in perturbation expansion as functions of the cor- 
relation lengths {^} = the set of quartic cou- 
plings {u} = {u±,ui\,Ux}, the set of asymmetric cou- 
plings {7} = {7x,7|j}, and the wave vector modulus k. 
The parameters in the order parameter vertex functions 

o (JV,L) 

rQi - ajv-ii - ii (foi' ll^c notation see Appendix A in pa- 
per I) via relations ([7]) - (O combine the corresponding 
parameters of the multicritical GLW-model ([5]). Thus 
all order parameter vertex functions calculated with 
have the property 



.,^(U},fc,W,{7}) 



(21) 



m,k,{u}). 
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meaning that they are identical to corresponding func- 
tions of the multicritical GLW-model ([5]). For this rea- 
son no distinction between the correlation lengthes en- 
tering the left and right hand side of (PT|) is necessary. 
The correlation lengths arc defined from the two-point 
order parameter vertex functions at the left side with 
©J and on the right side with ^ (see Eqs.(A7) and 
(A8) in paper I). Vertex functions of the secondary den- 
sity can be expressed as functions of {u} instead of 
{u} by using ([T])-®- Especially the two-point function 

Trnm = {"m-o '^o)^^! which will be of interest in the fol- 
lowing, can be written as 



fmm{{0,kAu},{l}) 



m'k,{u},m (22) 



B. Dynamical model 

The dynamical equations of model A in paper II have 
now to be extended by appending a diffusion equation 
for the secondary density. One obtains 



±0 
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(23) 
(24) 

(25) 



In addition to the two kinetic coefficients r_L and Fn of 
the order parameter in the corresponding subspaces, a 
kinetic coefficient A of diffusive type for the conserved 
secondary density is now present. The stochastic forces 



^</>_l; and 



fulfill Einstein relations 



{0^^{x,t) e^^^{x\t'))^2t^S{x-x')S{t-t')S"^ ,(26) 
{9^^^^{x,t) e'^y,t'))=2r^^6{x-x')S{t~t')S^^ , (27) 
{0m{x,t) 0m{x' ,t')) = -2XV^S{x ~ x')S{t ^t') ,(28) 



III. RENORMALIZATION 

A. Renormalization of the static parameters 

As a consequence of the discussion at the end of subsec- 
tion III Al all vertex functions will be expanded in powers 
of the quartic couplings {u} of the multicritical GLW- 
model and the asymmetric couplings {7}. The renor- 
malization scheme introduced in section III in paper I 
remains valid and will be used in the following. The cor- 
responding definitions and relations can be found therein 
and will not be repeated here. In particular, we imple- 
ment the minimal subtraction RG schemeiSiil directly 
at d = 3 to the two loop order. In the current extended 
model additional renormalizations for the secondary den- 
sity mo and the asymmetric couplings {7} have to be con- 
sidered. The renormalized counterparts of the secondary 
density and the asymmetric couplings are introduced as 



Too 

7 



= ^mTO 



1/2 



(30) 
(31) 



where k is the usual reference wave vector modulus and 
£ = A — d. The geometrical factor Ad and the diagonal 
matrix has been defined in Eq.(8) and Eq.(17) of 
paper I. With ([50]) and (pij) at hand, the renormalization 
for the CD-CD two-point vertex Tmm function readily 
follows 



^ mm — ^m mm i 



(32) 



The additional Z-factor Z^. and the matrix are re- 
lated to the known renormalization factors of the multi- 
critical GLW-model as a consequence of the reducibility 
of the extended model to the multicritical GLW-model. 

From the condition that (fT9|) is also valid for the renor- 
malized counterparts of the appearing quantities the re- 
lation 



Z^Zr 



Z^„ Z A, ■ Z A 



(33) 



with indices a, /3 = 1, . . . , n± and i, j = 1, . . . , n\\ corre- 
sponding to the two subspaces. The dynamical two-point 
vertex function of the secondary density has a general 
structure quite analogous to the corresponding functions 
of the order parameter (see Eqs.(6) and (7) in paper II). 
One can write 



mm 



(29) 



where Tmm{{^}:k) is the static two-point function dis- 
cussed in the previous subsection and (imm {{£,}, k, oj) is a 
genuine dynamical function^, m is the auxiliary density 
corresponding to m. For shortness we have dropped the 
couplings and kinetic coefficients in the argument lists of 



follows. For the second equality relation (18) of paper 
I has been used. The Z-factors of the asymmetric cou- 
plings are determined by the renormalizations of the sec- 
ondary density and the 0^-insertions in the multicritical 
GLW-model. 

Relation (llSp establishes a connection between the cor- 
relation functions of the CD and the i/i'^-insertions in the 
multicritical GLW-model. This relation should be invari- 
ant under renormalization. Thus the renormalization of 
the secondary density is related by 



Z„2 = 1 + 7^.A(M).7 



(34) 



to the additive renormalization A({m}) of the correla- 
tion function of the 0^-insertions p7p in the multicritical 
GLW-model introduced in Eq.(15) in paper I. 
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B. Renormalization of the dynamical parameters 

The general form of the rcnormahzation of the auxil- 
iary densities (j)±Q , (/)||o, and the kinetic coefficients T± 
and has been presented within model A in subsection 
III A in paper II. It remains valid and will be used in the 
following. Of course new contributions occur to the dy- 
namical renormalization factors especially of the kinetic 
coefficients 



(35) 



due to the asymmetric coupling 7. 

Within model C additional renormalizations only are 
necessary for the auxiliary density rfiQ and the kinetic 
coefficient A. Thus we introduce 



mo 



X^Z^\. 



(36) 



In the case of conserved densities the dynamical function 
^^mm({C}: (|29p does not contain new dimensional 

singularities. Therefore the corresponding auxiliary den- 
sity mo needs no independent renormalization. The Z- 
factor Zm is determined by the relation 



Zm, Z .„ 



(37) 



Due to the absence of mode coupling terms the renor- 
malization of the kinetic coefficient A is completely de- 
termined by the static renormalization and Z\ is 



Z\ — Z„ 



(38) 



are taken at fixed unrenormalized parameters. Inserting 
the two loop expression of ({u}) (see Eq.(31) in paper 
I) we obtain 

Cm(M,{7})-^7i + ^7^- (40) 

The column matrix of the /3-functions for the asymmetric 
coupling 7 is defined as 



/3;(M,{7})^«:^ 



(41) 



Inserting Eq. ([3T|) into the above definition one obtains 
together with relation (|33p the expression 



/37(W,{7})- [(-f + u)i + C0^(W) 
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(42) 



There 1 denotes the two dimensional unit matrix. The 
matrix C02({u}) has been introduced in paper I (see 
Eq.(22)). The (^-function C,m is exactly known from (|39|) . 
Thus finally we arrive at 



/37(M,{7}) = 



e 1^1 

h -7 

22' 



B0^(M)-7 1 



+Cl2(M) 



•7. (43) 



The above expression is valid in all orders of perturbation 
expansion. ({u}) and C02({u}) are calculated in loop 
expansion within the multicritical GLW-model. Their 
two loop expressions have been given in Eq.(31) and (23)- 
(26) in paper I. 



IV. C - AND f3 - FUNCTIONS 

As already mentioned in the preceding section the 
renormalization of the GLW-functional remains valid. 
This validates also all ^- and /3-functions introduced in 
section IV in paper I. We do not repeat them here, al- 
though they will be used in the following. 

A. Static functions 

Apart from the three /3-functions /3„^, /3u^ and Pu^, 
and the two C-matrices B^2 and ^^2 appearing in the 
multicritical GLW-model (see section IV in paper I), 
an additional ^-function and a column matrix of f3- 
functions for the asymmetric coupling (|16p have to be 
introduced. The relations between the renormalization 
factors discussed in subsection IIII Al give rise to corre- 
sponding relations between the and /3-functions. It 
follows immediately from ([34|) 



C™(M,{7}) 



d In Z„ 
din K 



1 



7^-B^2(M).7, (39) 



where B^2[{u}) has been defined in Eq.(30) in paper I. 
The K-derivatives, also in the following definitions, always 



B. Dynamical functions 



Using relation ([38]) the ^-function Ca corresponding to 
the kinetic coefficient A is simply given by 



Ca(M,{7}) 



dliiZ 



X 



din K 



2Cm({^},{7}) 



(44) 



The dynamical ^-functions of the kinetic coefficients of 
the order parameter are defined by 



4?(M,{7},W)^-Tr^ 

dluK 



(45) 



In the model C dynamics, they get non-trivial contribu- 
tions from the asymmetric couplings j± and 7||. They 
read now in two loop order 



n\\ w±7±7|| 
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4?(M,{7},M)) = C^'^"n"ih7ii,^ii) 




^(l + 2w)ln- 



l + 2v 



+ Cr'^^("ll'"x''f') 



(47) 



where we have defined the timescale ratios 

r± Til 

The ratio v is equally defined to paper II as the ratio 



(48) 



II 



w± 



(49) 



and is therefore a function of w± and wy . In (I46|) and 
(|47p several C-functions of known subsystems already has 

been introduced. Cp^'' (ua, Ux , with a = || or ± are the 
(^-functions of the full multicritical model A presented 
explicitly in paper II (see Eqs (14) and (15) therein), 



JA) _ n± + 2 
^ 36 
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2wln 



2{l + v) ^ {l + v) 
m ■ 



1 + 2t; 



1 + 2i; 



where (m, 7, w) are the genuine (^-functions of model 
C within the rta-component subspaces without pure 
fourth order coupling terms (pure model A terms). They 
have been given explicitly ii)^ for a n-component sys- 
tem in two loop order. These contributions of model C 
in the n-component subspaces without the corresponding 
model A terms are^ 



c 



(c„) 




2 / 4 
—uh - 3 In 



(l + 2w) jl + w)^ 
1 + w 1 + 2w 

(52) 

The /3-functions corresponding to the timescale ratios 
(l48)l and ((49l) can be expressed in terms of the corre- 
sponding C-functions of the kinetic coefficients: 



dv 
dn 



(53) 



/3iO|| 



dn 

dw± 



Ca) , 



(54) 
(55) 



with the K-derivatives taken at fixed unrenormalized pa- 
rameters. 

Note that these equations are not independent but one 
of the three equations can be eliminated by the relation 
v{l) = w^^{l)/'w±{l), which of course holds also for the 
initial conditions. 



V. FIXED POINTS AND THEIR STABILITY 
A. Static fixed points 

The FPs of the couplings Ua and their stability have 
been studied in section V of paper I. Their values and 
the corresponding transient exponents have been listed 
in Table I there. Let us recall that depending on the 
values of ny and n± one of the following FPs is stable and 
governs multicritical behavior: the isotropic Heisenberg 
FP H(n|| -I- n^) with u*^ = u*j_ = u^, the decoupHng FP 
V with u| ^ 0, 7^ and u* = 0, and the biconical 
FP with -it|| ^ 0, u*j_ ^ 0, and u* ^ 0. For each of the 
FPs of paper I one can now determine the FP values for 
7^ and 7|| from the Eq. 



/37(K},{7*}) 



e 1 
'2 + 2 



+Cl=(K}) 



■ 7* = 0. 



7*^.B^.(K})-r )1 



(56) 



This splits each FP of paper I into a set of FPs in the 
combined {ii}-{7}-space, which are equivalent in statics 
but different in dynamics. Given the formula for 
fij, one can see, that Eq. (|56p includes two equations 
which have to be solved for given FP values of the quar- 
tic couplings {u*}. The static FPs of paper I can be 
roughly separated into three classes: i) the Gaussian FP 
with u* = for all couplings; ii) the decoupling FPs 
7Y(n_L), 7Y(n||), 2? where = 0; iii) the isotropic Heisen- 
berg and biconical FPs 'H{n± + n||),i3 where all u* are 
different from zero. Henceforth we list the FPs values 
of the corresponding asymmetric couplings 7||' and 7^ in 
Tab. m that summarizes our analysis given below. Note, 
that 7^ = 0, 7^^ = is of course always a solution of 
equation ((56|) . independent which values {u*} have. We 
do not list this trivial solution explicitly in Table |T1 al- 
though this may be the stable FP for definite values of 
n±_ and ny in some cases. 
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1. Gaussian fixed point Q 

At this FP one has ~ and the two equations in 
56l) reduce to the condition 



(57) 



which is vahd in aU orders of perturbation expansion. 
The above equation defines a hne of FPs. 



2. Heisenberg and decoupling fixed points 7i(nx), 7i(n|| ), 2? 

At these FPs, where the cross couphng Ux vanishes, 
the matrix C,jp has the form 



^^(Ux =0) = (nil). . 



(58) 



The function c}'^^"'^ (ua) is the weh known C- function of 
the nQ-component isotropic system. Eq. (j56p reduces to 



(^-| + lr^.B,.(K})-r + cJ"^(-p)7ii 



where the matrix 5^,2 is of the form 



B^2{Uy^ = 0) = 











(59) 



(60) 



(61) 



The non trivial FP values for and 7|| resulting from 
Eqs.lUll) and are listed in Table U They are valid in 
all orders of perturbation expansion. 

We want to remark that in Table |T] only FP which 
exist for arbitrary order parameter component numbers 
are given. For special n- values a line of FPs exist where 
both asymmetric couplings 7q are different from zero. 
In the case = n\\ = n/2 one has u*^ ~ u^^ = u* 
and the ^-functions in ([50]) and ([SO)) are equal leading to 
lf + lf = {s-2CP^\u*))/B^-''\u*). 



3. Isotropic Heisenberg and biconical FPs Ti,{n± + ny), B 

At the isotropic Heisenberg FP T-l{n±^ + n\\) and the 
biconical FP B, where all couplings Ua are different from 
zero, it is more convenient to transform the matrix 
into its diagonal form with the transformation 



C+ 
c- 



= ■ C 



(62) 



introduced in section VLB in paper I. Inserting (|62|) into 
Eq. ([56|l leads to the transformed /3- function 



(63) 



with 



h -7 

2 2' 



■ 7 



C+ 
C- 



7± • 



(64) 



Here, the transformed asymmetric coupling column ma- 
trix is defined as 



7± 



p-1.7. 



Note, that the scalar quantity 

7^ • . 7 = 7I • B^,f 



■ 7± 



(65) 



(66) 



where B]^' = -3^2 P, is invariant under transforma- 
tion. Therefore in (|64p it is written in the untransformed 
form. At the FP Eq. ([63|l reduces to the condition 



(67) 



P^J{u*},{l*})^0 



because the determinant of matrix P does not vanish. 
Subsequently, Eq. ([M)) leads to two FP equations 



-£ + 7*^-Bl2 •7* + 2C-l7- =0 



(68) 
(69) 



In the above equations we have introduced the short hand 
notations B'^2 = B^2{{u*}) and = (±{{u*}). If both 
transformed asymmetric couplings 7^ and 7* are differ- 
ent from zero, the above two equations lead to the con- 
dition = C- • This condition is not valid if all quartic 
couplings u* are different from zero. Thus at least one 
of the two transformed asymmetric couplings, 7+ or 7_, 
has to be zero at the FP. The transformation matrix P 
has been presented in Eq.(64) in paper I. Expressed in 
terms of the C-functions it reads 



P = 



Pii P12 
P21 P. 



22 



/ 1 

[C^]l2 



[C2] \ 



(70) 

where [C02]jj are the elements of the matrix (for the 
two loop expressions see Eqs.(23)-(26) in paper I). 
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FP 


7l/7^ 


if 


7r 


G 


line of FPs (ETJ 


line of FPs (I571l 


line of FPs ([STI 










■2-e 




oo 




"II 





oo 


2 

—e 

















oo 


4''' 









V 











2(e-2C* ) 2(e-2C*) 

--(it!-)" 



[^^^]21 



TABLE I: Fixed points of the asymmetric couplings 7x and 7|| of the extended 0(?i||) © 0{ni_) model. The values of the 
isotropic Heisenberg FP T-l{n± + ny) and the Biconical FP B are valid in two loop order because Ea. (|75[) has been used. For 
all other FPs the expressions are valid in all orders of perturbation expansion. 



Let us now consider the two cases where one of the 
asymmetric couplings is nonzero. 

Case a: 7^ 0, 7* = 0: Taking into account Eq. ([65]) 
the condition for a vanishing 7* reads — P2iJ± + Piil^^ — 
0. Given the matrix elements (1701) it can be rewritten as 



(71) 



7i a - [c;^] 



22 



At finite 7* the bracket in Eq. ([55|) has to vanish, which 
results in the condition 



2c = • BI2 • 7* = £ - 2a 



a 
v 



(72) 



The last equality uses the definition of the asymptotic 
exponents derived in paper I (Eq. (90) there). 
Case b: 7* = 0, 7! 7^ 0: In this case Eq. ([S5)) leads im- 
mediately to the condition P227I — -P'i27j|' = 0- Inserting 
(|70)) gives 



71 



7f C-[0]n 



(73) 



At finite 7* the bracket in Eq. ([M|) has to vanish, which 
results in the condition 



2C = ■ BI2 • 7* 



2C* = 2- - d 
v 



(74) 



Again, the last equality uses the definition of the asymp- 
totic exponents derived in paper I (Eq. (82) there). The 
above equations ([7T|) . (|72|) and ((73)) . ((74)) respectively, de- 
termine the FP values of the two asymmetric couplings 
in the corresponding cases. The relations are valid in all 
orders of perturbation expansion. 

Since in two loop order B*^2 is diagonal and indepen- 
dent of the couplings {u} (see Eq.(31) in paper I) the left 
hand sides of Eqs ([72]) and read 



7 



*T 



2 • 7 



*2 ' 



^7± 



-7l| 



*2 



(75) 



In consequence the asymmetric static couplings are zero 
when the exponent expressions on the right hand side of 
([7^ and ([7^ are zero. This is the case if the specific 
heat like CD and/or the magnetic like CD susceptibility 
do not diverge. 

Using ([TS]) together with dill)-!!!!]) leads to the FP val- 
ues of the asymmetric couplings ^'^^ and 
Case a: 7* ^ 0, 7* = 0: 



if 



+ nil 



(76) 
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FIG. 1: Exponents ajv and 20/;/ — 3 appearing in Eqs. (I72p 
and (|74l) at ny = 1 as function of nx- Note that via relation 
(110511 2(f>/u — 2) = z — 2 the strong dynamical scaling exponent 
z results. 



FIG. 2; Fixed point values of the asymmetric static couplings 
7ll and 7x for Case b (7^ =Q, 'y*^ ^ 0) at the Heisenberg FP 
(nx < 1.61) and the biconical FP (1.61 < nx < 2.18). 



„*2 



2(e - 2C; ) 



(77) 



Case b: 7* = 0, 7I 7^ 0: 



if 



2(e-2C) 



(78) 



^11 



2(e-2C) 



[c; 



(79) 



Note that the ratios in ([71]) and ([73]) might be negative 
leading to a negative product 7j|7l- 

The explicit values of the above FPs depend on wether 
the isotropic Heisenberg or the biconical FP is inserted 
into the ^-functions. Eqs. (|76|) -(|79 l) are valid up to two 
loop order. In three loop order it is known from the 
isotropic GLW-model that the function B^2 gets v?- 
contributions"'^'^. In the multicritical GLW-model the ma- 
trix B^2 may also be non diagonal, and then Eg. ([75)) does 
not hold in this simple form. 

In the case of the isotropic Heisenberg FP u\ = 
Mjl = w* = u* Eqs. ([75|) - ([75|) simphfy considerably. The 
ratios of the elements of the ^,42 -matrix reduce to 



I 12 



7l 



(80) 



for Case a and 



2I 



(81) 



for Case b. For the second equalities ([TT]) and ([75)) have 
been used. Together with the relations e — 2(^;^ — ajv and 
£ — 2C* = — d, which introduce the critical exponents 
and follow from Eqs. (80), (81) and (82) in paper I, the 
values for the isotropic Heisenberg FP are: 
Case a: 7* 7^ 0, 7* = 0: 



(82) 



Case b: 7^ = Q, 7* 7^ Q: 



7f = 



T|| 



*2 



n± + nil n_L V ^ 



n± + n|j nil V 1^ 



(83) 



(84) 



Note that due to the sign in Eq. ([5T)) the relation 

, . 2 



7|i7± 



n± + ri|| \ I' 



2--d 



(85) 



holds in this case. For 7i|| = 1 and n± = 2 our results 
agree with those of Ref.^. 
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4- Resummation procedure 

As in our former papers^i^ of this series, in order to 
get numerical estimates we proceed within fixed dimen- 
sion RG technique, i.e. we evaluate RG expansions in 
couplings {m||,uj^,Ux} at fixed d = 3. Furthermore, as 
far as the expansions are known to have zero radius of 
convergence we use resummation technique^^ to get reli- 
able numerical estimates. The results given below were 
obtained within such a technique applied to the two-loop 
RG expansions. One of the ways to judge about typical 
numerical accuracy of our data, is to give an estimate 
for some cases where the expansions (and, subsequently, 
their numerical estimates) are known within much higher 
order of loops. As far as the static exponents a and v 
explicitly enter many of formulas considered above, let 
us take them as an example. Namely, let us estimate 
relations 



(20/z/-d)|d=3 

Ci/v\d=3 



2/i/_ 



(86) 
(87) 



that enter the formulas for the couplings 7^, 7||. The 
exponents and v- have been defined in Eqs. (80) 
and (81) of paper I. Fig. [T] shows the dependence of 
icji/u — 3 and of ajv on the order parameter component 
numbers n\_ at fixed rty = 1. Recall that of main in- 
terest for us will be the physical case n|| = 1, = 2 
indicated by the arrow. The region of shown in the 
figure covers also the region of stability of the Heisenberg 
0(n)-symmetrical FP, with n = n|| -Vn\^. In particular it 
starts at = 1 near the marginal field dimension ric at 
which the exponent a changes its sign. For the 0(n) vec- 
tor model an estimate based on the fixed d = 3 six loop 
RG expansions readai^: Uc = 1.945 ± 0.002. We get for 
the correlation length critical exponent in the Heisenberg 
0(2) FP: I' = 0.684, via hyperscaling relation this leads 
to a = —0.053. Our estimate correctly reproduces the 
absence of a divergency in the specific heat of 0{2) model 
(a is negative), however the value of nc ~ 1.6 we get is 
rather underestimated. Note however, that the fixed d 
approach we exploit in two loop approximation is essen- 
tially better than the corresponding e expansion. Indeed, 
in two-loop e-expansion one gets: Uc = 4 — 4e, which does 
not lead to reasonable estimatesi^. The two loop esti- 
mate of the massive field theory at d = 3, nic — 2.01—, 
is more close to the most accurate value of Ref.^-, how- 
ever it gives a wrong sign for the exponent a. In any 
case the negative value of a for n±_ = 2 agrees with other 
calculations as reported in paper I. 

It turns out (see below) that Case b is the stable FP 
for the asymmetric couplings. In order to evaluate nu- 
merically the values of couplings 7^, 7^^, we therefore 
substitute the resummed fixed point values of the static 
couplings ,u±,Ux} into formulas ([75)) . I|79p and resum 
the resulting expression. In principle, one can use differ- 
ent ways for such an evaluation. Indeed, as we proceeded 
before, one can present these formulas in the form of ex- 
pansion in renormalized couplings (keeping the two-loop 



terms) and resum the resulting second-order polynomial. 
Alternatively, based on the observation that numerator 
and denominator of Eqs. ([78|) . ([79|) contain combinations 
of critical exponents, one can resum the numerator and 
the denominator separately. We will exploit both ways 
which naturally will lead to slightly different numerical 
estimates. This difference may also serve to get an idea 
about typical numerical accuracy of the results. Sepa- 
rately, we will evaluate the ratios 7±/7||. Again, it will 
be done by resummation of the series for this ratio, Eq. 
(|73l) . as well as by using resummed values for 7j|, ^'j_.ln 
particular, for n|| = 1, = 2 we get: (7^)^ = 0.034, 
(7||7kr)^ — 0.286 (when denominator and numerator are re- 
summed separately), and (7^)^ = 0.031, (7||*)^ = 0.293 
(when an entire expression is resummed). Resulting dif- 
ferences of the order of several percents bring about a 
typical numerical accuracy of the estimates. In Fig. [2] 
we plot the FP values of the asymmetric couplings and 
their ratio obtained within resummation of the entire ex- 
pressions. These values will be used below to calculate 
the critical dynamics. 



B. Static transient exponents 

The stability of the fixed points is determined by the 
sign of the corresponding transient exponents. Latter can 
be found from the eigenvalues of the matrix 



97/3 



with a,f3 =_L, ||. Inserting (|43p into (|55|) the correspond- 
ing eigenvalues read 



± 



'7i_L7j - n||7f? 



22 



-,l/2> 



+ (ji^7j.7|| +2[g2]^2)(n||7j.7| +2[C^2]2i) 



.(89) 



The above eigenvalues are valid in two loop order because 
Eq. ([75|) already has been used. The transient exponents 



;(±)^A(±)(m = K},{7} = {7*}) 



(90) 



are calculated by inserting the fixed point values of the 
static couplings into the eigenvalues. With the fixed 
point values ([82]) - and Eq. ([Ml) we obtain for the 
isotropic Heisenberg FP the transient exponents 
Case a: 7* =^ 0, 7* = 0: 



(+) 



(-) 



(91) 
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Case b: 7^ = 0, 7I 7^ 0: 

^(+)^2^-d, (92) 

V 

W* is the root 

^=\/([C0^]ll-[<0^]22)'+4[C0^]l2M21 (93) 

taken at the fixed point values of the couphngs. It is 
always positive and reads at the isotropic Heisenberg FP 
in two loop order 

= (1 - y) • (94) 

Thus one concludes that Case b is the stable FP even if 
a would be positive^^. For the biconical FP the stability 
of Case b can be verified explicitly by the flow of the 
couplings. 



C. Dynamical fixed points 

Calculations of the dynamical FPs values are done by 
solving the FP equations for the dynamical /3-functions. 
Since only two of the equations ((53)) - (l55|) are independent 
the third equation serves as consistency check of the so- 
lution found. It is useful to choose for this purpose Eqs 
([54)1 and ([55]) for the time scale ratios and w±. Then 
one has to solve 

/3u,ii(w||,wj^,W||/w±) = 0, (95) 
/3„^(w||, wx, W||/w±) = 0. 

To find the dynamical FP values coordinates, the re- 
summed FP values of static couplings U|J , , u* , 7^ , 7|* 
are inserted into these equationsi^. 

The dynamical FPs depend on which static FP is con- 
sidered. There might be several dynamical FPs for one 
static FP, which could be either strong dynamical or weak 
dynamical scaling FPs. Since also unstable static FPs 
might be reached in the asymptotics if one starts with 
static initial conditions in the attraction region of this 
FP (a subspace in the space the static couplings see e.g. 
Fig. 3 in paper I) at least both the Heisenberg FP and 
the biconical FP have to be taken into consideration. It 
turns out that in both cases apart from the trivial unsta- 
ble FP where all timescale ratios are zero (see below) two 
dynamical FPs are found: (i) an unstable weak dynami- 
cal scaling FP corresponding to model A and (ii) a stable 
new strong dynamical scaling FP. In the physical inter- 
esting case n|| = 1 and n± = 2 these cases correspond to 
dynamical behavior at a multicritical point of bicritical 
and tetracritical type respectively. The different types of 
weak and strong dynamical FPs are shown in Tab. [Til 
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' 1 ' 1 ' 1 

Heisenberg FP — ^ 


biconical FP— •r 


1 


^ — strong dynamic 


1 
1 

1 
1 
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/ 
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D I . I ■ , 1 , 1 , 1_ 
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perp 

FIG. 3: Fixed point values of the timescale ratios v, l/w± 
and «)|| for the static stable FPs: the isotropic Heisenberg FP 
(nx < 1.61) and the biconical FP (1.61 < n± < 2.18). Strong 
dynamical scaling is valid up to the stability borderline to the 
decoupling FP. In the biconical region the values of vc and 
1/wcx are finite but not to distinguish from zero on this scale. 
The notation of the dynamical FPs correspond to the notation 
m Tab[ni The dashed curve shows the unstable model A FP 
(see text). 



1. Strong dynamical scaling fixed point 

In order to find the strong dynamical scaling FPs it is 
not necessary to discriminate between the static Heisen- 
berg FP or the biconical FP although the dynamical 
equations to be solved simplify in the first case a little 
bit. Thus we use the results for the FP values derived in 
paper I for the quartic couplings {u} and the FP values 
for the asymmetric couplings {7} of Case b (see Eqs. 
([73| . ([78)1 and ([79)) ). At the strong scaling dynamical 
FP all timescale ratios have to be nonzero and finite. 
Moreover due to the definitions of the timescale ratios it 
follows that W|| = v*w\. This dynamical strong scaling 
FP value is found by setting the differences of two of the 
three dynamical C-functions, ([44 |) - ([47|) to zero leading to 
three equations. Since there are only two independent 
time ratios the third equation can be used to check the 
results. 

The FP values (FP with subscript C in Tab. [IJ) have 
been plotted in Fig. [3] for different n± at n|| = 1 and 
the numerical values for n±_ = 2 are collected in Tab. 
mil This shows that the FP value of the timescale ratio 
for V is different from the FP value found in the pure 
relaxational model A. These were at the Heisenberg FP 
v\ = 1 and at the biconical FP v* = with — > 00 
in approaching the stability borderline to the decoupling 
fixed point (see Fig. 1 in paper II). 

A numerical problem arises in finding the FP values of 
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FP 



scaling type 





weak 







oo 


Ha 


weak 










He 


strong 









2^ - 1 
2 + 0] 
2i- - 1 



infinite fast 
2 + CJ7 
2^ - 1 



2± 
2± 



1 
1 

2^ - 1 



Be 



weak 
weak 
strong 



(B) 
„(B) 







CX) 





2^^ - 1 
2^ - 1 



infinite fast 
2^ - 1 



2^ 
2^ 



TABLE 11: Types of dynamical FPs for the static Heisenberg and biconical FP. Not included is the trivial unstable fixed point 
with all times scale ratios equal to zero. The value of c reads c = 61n(4/3) — 1. Note that for the weak scaling FP the result 
is only valid in two loop order, whereas the relation of the dynamical critical exponent in the strong scaling FP holds in all 
orders. 



FP 










if 


V* 


™1 


wl 


Be, a 


1.28745 


1.12769 


0.30129 


0.29378 


0.03170 


6.09592 • 10'*^ 


1.24285 - 10*^ 


0.75763 


-He, b 












7.29393-10-^ 


1.55665-10'' 


1.13541 


1.00156 


1.00156 


1.00156 


0.72554 


0.18139 






1.13541 


He, a 


7.30771-10-^ 


1.55372-10'' 



TABLE III: FP values of couplings and timescale ratios for n\\ = 1, n± = 2. a FP values of the timescale ratios found via 
approximation using Eqs. (100), (101), as described in the text with the values for FP B: A = 0.09770, B — 0.00101 and the 
FP Ti: A = 0.31534, B = 0.03311. b numerical solution for the FP values of the timescale ratios; 



V and 1/wj^when they reach very small values. It cannot 
be numerically decided wether the FP values are zero or 
finite. In order to clarify the existence or nonexistence of 
a weak scaling FP one has to look for an analytic expres- 
sion for the small FP values. However it is numerically 
easy to find the FP value of W|| which is nonzero and finite 
in the whole region up to the stability borderline between 
the biconical and decoupling FP. In order to solve this 
problem the dependence of the (^-functions are studied 
within this region. One observes that there are logarith- 
mic terms which would diverge in the limit u — > under 
the condition w±^ = w\\/v. Thus one obtains two equa- 
tions for the FP value of v and w\\. In the equation for 
the FP of one might safely perform the limit v ^ 
and w \ oo. This leads to 







4?HK},{7*},^ = 0,u;||, 



oo 



2C- (96) 



Using limiting functions 



+ 2 ^„ / 4 1 



72 



(97) 



and 



4 



*2 



1 -I- U^ll 



1 - - 



4( 



1 - 3 In - 



3 



W||7f / nil 



(98) 



3(r 



11+2)^^4 {l + 2w\\) ^^ {l + w\\Y 
2 ^3 l + wii ^ H-2wii 



then Eq. (IM]) reads 




2 / 4 
—uti 1 - 31n- 

II V 3 



W\\ 



3(n| 



+ 2) 4 

' m - 

2 3 



l + 2u;|, (1 + W||)2 
m - 



1 + 2w 



(99) 



^117117^ 



36 



^11 



4 1 



3 In - - 



f"||7||7± 

1 + W\\ 



72 



2C 



This equation is solved numerically to give the value of 
uiji which then is inserted into the second equation for v* . 
In order to find v* one collects the logarithmic diverging 
terms in the equation for v* , 



(^f,\'",'^\\iW\\/v) - C'^^{v,w\\,w\\/v) 



iC), 



0. 



(100) 



In the remaining terms the limit w — > can safely per- 
formed. Then the solution reads 



In v* 



A 
B 



(101) 
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FIG. 4: Dynamical flow at — 1 and different n± values for 
different dynamical initial conditions numbered 1 to 4. The 
static couplings are chosen to be fixed at their stable FP values 
(the isotropic Heisenberg FP for n± — 1.2, the biconical FP 
for n± — 1.7). The dynamical FP values are u* = 0.399, 
0.004, W|| = 1.661, 1.300 and wl = 4.159, 351.06 at He, Be 
respectively. Also shown is the surface v = •u}\\lwx_ to which 
the flow is restricted. 



with 




n_L + 2 
3 



<(l 



31n- 



FIG. 5: Dynamical flow at n\\ — 1 and nx ~ 2 for different 
dynamical initial conditions numbered 1 to 4. The static cou- 
plings are chosen to be fixed at their biconical FP values. The 
static and dynamical FP values of Be are given in Tab. IIIII 
The dynamical FP lies outside the region shown. Also shown 
is the surface v — w^^/w± to which the flow is restricted. 



\ 2 2 3 2 



7^711 



In 2 u I 

36 ^ 



and 



/ ''ill o n\i 



(103) 



It can be shown that A and B are positive. Approach- 
ing the stability borderline to the decoupling FP A stays 
finite and B goes to zero since and 7^ go to zero. 
In consequence v* goes to zero and w* goes to infinity 
in and only in this limit. The analytic solution found 
within this region joins smoothly to the numerical solu- 
tion found for larger values of the timescale ratios. Thus 
it is proven that in the whole region where the Heisen- 
berg FP or the biconical FP is stable dynamical strong 
scaling holds. 

Considering the FP values for the timescale ratios for 
the Heisenberg FP in the region of > 1.7 (where it is 
reached only for static initial conditions in a subspace of 
the fourth order couplings) one also finds a small value 
for V* , a very large value for w* and a nonzero finite 
value for w*^. However contrary to the biconical FP now 
A and B stay finite at the stability borderline between 
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the biconical FP and the decouphng FP at n± ~ 2.18. 
Indeed the values calculated for the Heisenberg FP from 
Eqs (IMl) and are A = 0.34866, B = 0.02558 

whereas for the biconical FP one obtains A = 0.05268, 
5 = 4.27958-10-9. 

The asymptotic dynamical exponents are obtained 
from the values of the C-functions at the FP: 



z^ii — 2 + Cr|, 



2 + 



2 + Ca. (104) 



At the strong dynamical scaling FP all dynamical 
functions are equal to twice the static C- function (^m- 
Therefore the CD induces the value of the dynamical 
critical exponent z 



2C™ 



1 



(105) 



according to Eq. ([39l) and Eq. ([74]). The values for the 
static exponents depend on which static FP is stable. For 
nil = 1 the n^-dependence of z — 2 is shown in Fig. [1] 



2. Weak dynamical scaling fixed point 

Weak dynamical scaling FPs are solutions of the dy- 
namical FP equations where one or more of the FP values 
of the timescale ratios are zero or infinite. Such a weak 
dynamical scaling FP has already been found in model 
A and it became stable at the stability borderline to the 
decoupling FP. 

Indeed Eqs ([95]) allow solutions where both timescale 
ratios w\\ and w± are zero. In such a case one has to rely 
on the third equation for the ratio v = w^^/w^ to find 
the limiting FP value. However in the limit W|| — > and 
w± — > Eq. ((53|) for v reduces to the FP equation of 
model A (FP with subscript A in Tab. [n|. Thus one 
recovers the model A FPs in this case. 

There is no solution w^^ = v* — and wj_ nonzero and 
finite due to the Inv term in (jSOp . For a similar reason 
no FP with w*^ nonzero and finite, — and v* = oo 
is possible. However a FP with wJl nonzero and finite, 
Tj* = and w± = oo is possible (FP with subscript in 
Tab. |TT|). The values of w| are obtained from Eq. ([M)) . 
but now these values are not an approximation but the 
exact Ciu FP values for any n\\ and n±. 

The dynamical critical exponents may be different in 
the case of weak dynamical scaling. For the weak model 
C FP (subscript C^) w*,^ is finite and nonzero therefore 

= Ca and Z0 = 2^ - 1. Thus the CD sets the 



timescale for the OP 



'II- 



Inserting w*^_ ~ oo into Cr^ 



leads due to logarithmic diverging terms to an infinite 
value of the corresponding dynamical exponent z^j_ . This 
indicates that the density (j)± is much faster than the 
other densities. It is especially much faster than the other 
OP 0)1. 

In the case where both FP values of the timescale ratios 
wn and wj_ are zero and v is finite and nonzero, both OPs 



have the same timescale with the dynamical exponent (of 
model A) z = 2 + crj different from the exponent of the 
CD Zm. =2^-1. 



D. Dynamical transient exponents 

The dynamical transient exponents can be calculated 
from the matrix of the derivatives of the /3-function with 
respect to the timescale ratios u, w\\ and wj_. Since only 
two timescale ratios are independent only two are con- 
sidered in the stability matrix. The eigenvalues of the 
2 X 2-matrix have to be positive for the overall stable 
FP otherwise the FP is unstable. In the following the 
timescale ratios and w± are chosen as independent. 

The model A type FP with v* = va nonzero and finite 
is unstable since the two eigenvalues. 



— Cr,, 2C*j 



and 



=C?, -2C (106) 



are negative. In fact they are equal because Cr = Cr^ ■ 
Their values are calculated with 2(*^ — 2(j)/v ~ 1 and 
inserting the model A FP value v* = va (see Fig [3]) and 
w| = = 0. ThenC?||-2(/)/:/-Hl < Oandw^n = -0.126 
for n|| = 1 and ny = 2. 

Similarly the instability of the FP with v* ~ wj| = 
w* = can be shown. However some care has to be taken 
due to the vanishing timescale ratio v. The eigenvalues 
are again given by Eq. (|106p but now they are different. 
Whereas lOw^^ is negative goes to oo due to the Inw 
term in the model A C, function ([50)1 . 

The transient exponent for the strong scaling FP are 
the eigenvalues of the matrix of derivatives of the /3- 
functions according to the timescale ratios at the FP 





dww 








dww 






(107) 

Use has been made from the independence of (\ on the 
timescale ratios. The nondiagonal elements depend on 
the timescale ratios by which they are derived only via 
V and therefore are proportional to l/w±. In the region 
where w± is very large the two eigenvalues are then given 
by the diagonal elements. 



dw 



... fdCr 



(108) 



Near the stability borderline to the decoupling FP the 
second eigenvalue goes to zero according to 



B + 0{l/wA_) 



(109) 



with B from Eq. (|103p . being exactly zero at the bor- 
derline. The value of the slow transient at ny = 1 and 
nj^ = 2 is given by = 0.001 Thus as shown in the 
next section in this case nonasymptotic effects are present 
in the physical accessible region. 
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As already mentioned for the Heisenberg FP Tic B 
does not reach zero at the borderUne but its value is 
one order smaller than the static transient exponents. 
At n±^ — 2 and n±^ = 2.18 the value of the dynamical 
transient exponent is given by to^^ = 0.033 and uJw^ = 
0.026 respectively. 



VI. DYNAMICAL FLOWS AND EFFECTIVE 
EXPONENTS 

The flow of the timescale ratios is described by the RG 
equations 



dv 



dl 



dw 



l- 



dl 



/3.(M,{7},M), 

PwA{u},{l]-,{w]) ■, 
PwA{u},{l},{w}) , 



(110) 



with the /3-functions Eqs. (153l)-(|54l)- Note, that the dy- 
namical /3-functions depend also on the RG equations of 
the static quartic couplings (Eqs (33)-(36) of paper I) 
and the RG equations of the asymmetric couplings 



;|I = ^,(M,{7}) 



(111) 



with the /3-function Eq. (l43l) . 

In order to simplify the picture it is assumed that the 
static couplings have already reached the FP by which 
they are attracted from their initial conditions. Then the 
RG flows are displayed in the three-dimensional space of 
the time-scale ratios w\\,w_l,v in Figs. S] and [5] for the 
physical interesting case at n\\ ^ 1^ — 2. The static 
biconical FP is stable for this case in general. However 
for initial conditions on the surface separating the stable 
biconical FP from the Heisenberg FP the flow is attracted 
to the Heisenberg FP. Therefore one may also fix the 
static parameters to this FP. 

To give an overview of the different patterns of the 
flows three different value of n_L are chosen for fixed n\\ — 
1: (i) = 1.2, where the static Heisenberg FP is stable 
(Fig. HI), (ii) nj. = 1.7 (Fig. m) and (iii) n^^2 (Fig. O 
where the biconical FP is stable. 

In all cases the FP values of the timescales are nonzero 
and finite but the values of v* becomes very small and 
■w\ very large. The asymptotic approach to the FP in 
cases (ii) and (iii) occurs in the direction of the w^-axis 
almost at z; and W| ^ w*y 



A. Effective exponents 

We define the effective exponents by: 

z|| ,//(/) = 2 + Cr„(M0},{7(0},M0}), 
z^e//(0 = 2 + Cr,(M0},{7(0},M0}), 

Zmeff{l) = 2 + U{u{l)}). 




2.00 
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FIG. 6: Effective dynamical exponents z\\, z±, and Zm calcu- 
lated along the different RG flows of Fig. [5] (indicated by the 
numbers). The insert shows that even for flow parameters 
as small as InZ = —2000 the effective exponent 2x has not 
reached its asymptotic value 2.18. 



These exponents appear e.g. in the critical temperature 
and/or wave vector dependence of the transport coeffi- 
cients describing the relaxation of the alternating magne- 
tization or the diffusion of the magnetization in the direc- 
tion of the external magnetic field. They are in principle 
experimentally accessible. An interesting feature is the 
independence of the effective dynamical scaling exponent 
Zm of the CD from the dynamical timescales. Therefore 
its nonasymptotic value is only due to nonasymptotic ef- 
fects within statics. This allows to trace back nonasymp- 
totic effects in dynamical quantities to the slow transients 
in statics or those appearing in dynamics. 

In order to calculate the numerical values of the ef- 
fective exponents, we substitute into Eq. (I112p the re- 



summed coordinates of the static FP {uZ 



, } and 



the values of timescale ratios (£), i(;j^(^), w(^)} along 
the RG flow. Choosing the FP values of the static cou- 
plings fixes the asymptotic values of the effective dy- 
namical exponents since they are expressed by the static 
asymptotic exponents for the strong dynamical scaling 
FP. 

The results shown in Fig. [S] correspond to the case 
n\\ = 1, = 2. We evaluate the timescale ratios along 
the previously obtained flows 1-4 in Fig. [5] The effective 
exponents for the different initial conditions are shown 
by numbered solid lines, As one can observe from this 
figure, the exponents calculated along several flows do 
not coincide for the values of the flow parameter shown. 
However one sees the merging of the different values for 



to their asymptotic value z*. 



2.18 given by the 



Jll2) static value corresponding to the CD (the constant line 
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in Fig. [6]). More remarkable is the difference between 
the effective exponents for the paraUel and perpendicu- 
lar components of the OP. This difference can be traced 
back to the fact that 'w± did not reach its very large FP 
value due to the slow dynamical transient. Only for flow 
parameter values larger than \nl ^ —2000 the effective 
exponent attains its asymptotic value z± ~ 2.18 as it 
should be (see the insert of Fig. [S] where the smallest 
value is In / = -2000 for the flow 1). 

Furthermore, we see that the exponent Z|| flows to- 
wards its asymptotic value Z|| = 2.18. In fact, exponent 
z± attains the asymptotic value z± — 2.18 as well, but 
for much smaller values of the flow parameter. In the in- 
sert of Fig. [6] we show this exponent for all flows within 
larger range of Ini*. 

For certain initial values of the static parameters it is 
possible that the unstable Heisenberg FP is reached (see 
Fig. 3 in paper I the flow number 1, which lies on the 
surface separating the attraction region of the biconical 
FP from the flow away solutions. In such a case the 
effective exponents reach faster the asymptotic value of 
the dynamical exponent z = 2.44 since the dynamical 
transient exponent lOw^ = B is lager by a factor of 33 
(for the values of B see caption to Tab. IIIip . 

VII. CONCLUSION AND OUTLOOK 

The effect of coupling the conserved magnetization 
parallel to the external magnetic field to the two OPs 



(the components of the alternating magnetization paral- 
lel and perpendicular to the external magnetic field) leads 
to strong scaling dynamical behavior. This means that 
the timescales of all dynamical quantities scale with the 
same dynamical critical exponent z — 2(j)/v — 2. However 
this scaling behavior might be hidden by nonasymtotic ef- 
fects dominating the physical accessible region when ap- 
proaching the multicritical point due to a very small dy- 
namical transient. This applies in the physical interesting 
case n|| = 1 and n± = 2 where one of the timescales is 
almost zero and the other one almost infinite. In con- 
sequence the magnetic transport coefficients might show 
different effective behavior with temperature when ap- 
proaching the multicritical point. The dynamical ampli- 
tude ratios might be far from their asymptotic values and 
show nonuniversal behavior. 



For a complete description in the whole n||-n_i_-space 
the model presented here has to be extended in two ways. 
First for the physical case n\\ = 1 and n± = 2 one has to 
introduce reversible terms in the equations of motions. 
Second one has to allow for an asymmetric coupling to a 
energy like CD in addition to the magnetization. 
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